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ABSTRACT 

Multiple or extended turnoffs in young clusters in the Magellanic Clouds have recently 
received large attention. A number of studies have shown that they may be interpreted 
as the result of a significant age spread (several 10 8 yr in clusters aged 1 2 Gyr), while 
others attribute them to a spread in stellar rotation. We focus on the cluster NGC 1856, 
showing a splitting in the upper part of the main sequence, well visible in the color 
niF 336 W- m F 555 W; and a very wide turnoff region. Using population synthesis available 
from the Geneva stellar models, we show that the cluster data can be interpreted as 
superposition of two main populations having the same age (~350 Myr), composed for 
2/3 of very rapidly rotating stars, defining the upper turnoff region and the redder main 
sequence, and for 1/3 of slowly/non-rotating stars. Since rapid rotation is a common 
property of the B-A type stars, the main question raised by this model concerns 
the origin of the slowly/non-rotating component. Binary synchronization is a possible 
process behind the slowly/non-rotating population; in this case, many slowly/non¬ 
rotating stars should still be part of binary systems with orbital periods in the range 
from 4 to 500 days. Such periods imply that Roche lobe overflow occurs, during the 
evolution of the primary off the main sequence, so most primaries may not be able to 
ignite core helium burning, consistently why the lack of a red clump progeny of the 
slowly rotating population. 

Key words: stars: early-type - (galaxies:) Magellanic Clouds - (stars:) Hertzsprung- 
Russell and colour-magnitude diagrams - stars: interiors - (Galaxy:) globular clusters: 
general 


1 INTRODUCTION 

Extended main sequence turnoff regions (eMSTO) appear 
to be a typical feature of massive intermediate-age glob¬ 
ular clusters (GCs) in both Magellanic Clouds, as shown 
by deep observations of their color-magnitude diagrams 
(CMDs) taken with the ACS and WFC3 onboard the Hub¬ 
ble Space Telescope (HST). Most of these clusters are in 
the age range 1-2 Gyr (Mackey et al. 2008; Milone et al. 
2009), and, if the eMSTOs are interpreted as due to dif¬ 
ferences in the formation epoch of stars, their age spreads 
range from 150 to ~500 Myr (Milone et al. 2009; Girardi 
et al. 2011; Goudfrooij et al. 2011; Rubele et al. 2013) . 
This possible age spread has been sometimes used (Con- 
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roy & Spergel 2011) as a possible evidence of the age dif¬ 
ference underlying multiple generations in old galactic GCs 
(Gratton et al. 2012; Piotto et al. 2012), although theoret¬ 
ical work on these latter shows that the possible age differ¬ 
ences must be contained within about 100 Myr (D’Ercole 
et al. 2010). Moreover, (limited) spectroscopic observations 
of stars in the eMSTO cluster NGC 1806 seem to suggest 
that this cluster has homogeneous abundance of O, Na, Al 
and Mg, and that the O-Na anti correlation may be not 
present in clusters with extended turnoffs (Mucciarelli et al. 
2014). If confirmed by larger dataset of spectroscopic studies 
in clusters having eMSTO, the result would indicate the lack 
of any connection between eMSTO clusters and multiple- 
populations GCs, since the anticorrelation between oxygen 
and sodium abundances is one of the key fingerprints of mul¬ 
tiple populations in old GCs. Notice that, in galactic GCs, 


© 2015 The Authors 


2 


F. D’Antona et al. 




Figure 1. Left: 300 Myr isochrones for Z=0.006 from Georgy et al. (2013) for non rotating models (grey line) and models rotating at 
cj= 0.5 (Jcrit (blue line); center: luminosity versus mass for the same isochrones; right side: time evolution of the 3 Mq. non-rotating and 
rotating models. 


the Na-O anticorrelation is present also at metallicities as 
large as those of LMC clusters, in which oxygen is slightly 
depleted, but the sodium abundance spans >0.5 dex (Car- 
retta et al. 2009). Other LMC intermediate age clusters may 
have sodium spreads by 0.1-0.35 dex, which may hint to a 
possible link with the presence of second generation stars 
(see the discussion in Goudfrooij et al. 2014), but the data 
(Mucciarclli et al. 2008) are as limited as those in NGC 1806, 
to allow final conclusions. A lack of the Na-O anticorrelation 
would call into question any attempt to use observations of 
the eMSTO young and intermediate-age clusters, to rule out 
formation and self-enrichment models for old GCs. 

Bastian & de Mink (2009) suggested that stellar rota¬ 
tion may cause the eMSTO, as stellar rotation affects the 
structure of the star and the inclination angle of the star 
relative to the observer will change the effective temper¬ 
ature, hence observed colour. This result has been subse¬ 
quently questioned (e.g. Girardi et al. 2011) with a num¬ 
ber of studies supporting the idea that an age spread is re¬ 
sponsible for the observed CMD features (Goudfrooij et al. 
2011; Rubele et al. 2013). Yang et al. (2013) emphasized 
that the evolution of rotating models drastically depends on 
the parameters adopted in the description of the rotational 
mixing, so that the eMSTO can be explained -or not- de¬ 
pending on which set of models are adopted. More recently 
Bastian & Niederhofer (2015) and Li et al. (2014) excluded 
age spreads in intermediate age LMC clusters, by inspec¬ 
tion of their sub giant branch morphology, but their analysis 
was challenged by a very careful study by Goudfrooij et al. 
(2015). At the same time, Brandt & Lfuang (2015) studied 
the turnoff spread expected on the basis of a rotation disper¬ 
sion in the component stars based on the Geneva database 
created by C. Georgy and S. Ekstrom 1 (Georgy et al. 2014) 
and showed that, in the color magnitude diagram of the 
typical HST bands mF435w versus color niF435w-mF8i4w, 
the turnoff area covered by turnoff stars increases, reaching 


1 http://obswww.unige.ch/Recherche/evoldb/index/ 


a maximum towards 1-1.5 Gyr, which are in fact the typical 
ages of clusters with eMSTO. 

In this paper, we focus our attention of the LMC cluster 
NGC 1856, which is much younger than any other cluster 
with evident eMSTO. The eMSTO has been recently iden¬ 
tified by Correnti et al. (2015) and Milone et al. (2015) who 
showed that this feature is consistent with an age spread of 
~80-150 Myr, although they do not exclude the rotational 
interpretation. 

Interestingly, in addition to the eMSTO, this cluster ex¬ 
hibits a split upper main sequence which is clearly visible 
when the F336W filter is used (like in the hifsssw versus 
mF 336 W-mF 555 w CMD). In a preliminary attempt dealing 
with models with age spread, Milone et al. (2015) show that 
this feature may be explained by models in which the stars 
belonging to the cooler main sequence have a larger metal- 
licity. While helium and light element variations are respon¬ 
sible for the MS splitting observed in the old galactic GCs, 
their effect on the double MS of NGC 1856 has been not yet 
investigated. In this paper, by using the Geneva database, 
we show that the main cluster features can be interpreted 
in the framework of the rotational hypothesis, and provide 
a possible explanation for the split MS. 


2 MODELS 

The rotational interpretation of the extended turnoffs is 
linked to the initial distribution of the cluster stars’ rota¬ 
tional velocities. Taking as an example the nominal age of 
300 Myr attributed to NGC 1856 by Milone et al. (2015), 
the turnoff stars have mass of about 3 Mq, so the hydrogen 
core burning phase is convective. The evolution of stars hav¬ 
ing a convective core in the main sequence is dramatically 
affected by rotation. In fact, the mechanism of chemical mix¬ 
ing associated with the transfer of angular momentum from 
the (more rapidly rotating) convective core to the stellar 
envelope provides fresh hydrogen to the burning core, and 
extends the main sequence phase. Fig. 1 illustrates the dif¬ 
ference in the isochrone location for an age of 300 Myr and 
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metallicity Z=0.006, having no rotation or a rotation equal 
to 0.5 the critical rotation rate (Georgy et al. 2013). The 
turnoff is more luminous by ~0.3 dex. The masses in evo¬ 
lution at the turnoff are ~2.9 Mq in the non rotating case, 
and ~ 3.2Mq for model rotating 0.5 cu C rit (central panel of 
Fig. 1). This example clearly illustrates how a range of rota¬ 
tional velocities in a population of coeval stars can be charac¬ 
terized by a broad range of turnoff luminosities and evolving 
masses, similar to those one would find in a cluster with non 
rotating stars spread over a range of ages. The hydrogen 
consumption and the luminosity evolution versus time are 
shown in the right panel of Fig. 1, for the 3Mq evolution. 
This panel makes clear that the time-luminosity evolution 
will depend on the detailed assumptions made concerning 
the mechanisms of chemical mixing. 2 

It is therefore important to keep in mind that the 
detailed relation between turnoff spreads and rotation rates 
needed to model the observed CMD features depends on the 
specific parameters adopted in the stellar models used; a 
smaller rotation rate range would be needed if, for example, 
a higher chemical mixing efficiency than adopted in these 
specific models were assumed. In particular, in the Geneva 
database we use, the turnoff luminosity is not linearly 
related with the rotation velocity assumed for the stars, 
while the effect of rotation of the MS T e fr is. We show this 
in Fig. 2, where isochrones of age t=350Myr are compared 
for oj /cUcrit-—0, 0.6, 0.8 and 0.9 in the theoretical plane. We 
see that a significantly cooler MS requires adopting a very 
large rotation. This effect is increased in the observational 
plane by the corrections due to limb and gravity darkening. 


3 THE COLOR MAGNITUDE DIAGRAM OF 
NGC 1856: COMPARISON WITH 
ROTATING MODELS 

Fig. 3 shows the CMD of the cluster NGC 1856 in niFsssw 
versus the niF 336 w-niF 555 w color observed by Milone et al. 
(2015) (left panel) and mF438W versus mF43sw-mF8i4w color 
(right panel). Fig. 4 shows the Hess diagram in the mF 336 W- 
niF 555 w color. This is the first evidence of an extended main 
sequence turn-off ( lmag) in a young cluster in the Mag¬ 
ellanic Clouds as discussed by Milone et al. (2015) and Cor- 
renti et al. (2015). 

In Fig. 4 note also the split of the upper MS in two com¬ 
ponent, which host ~ 33% and 67% of MS stars in the red 
and blue part respectively. The detailed analysis of Milone 
et al. (2015) has demonstrated that both the split of MS and 


2 These will be due both to the assumptions on rotational mixing, 
and to the form of overshooting of the convective core, adopted 
independently from rotation (that is, valid also in non rotating 
models). In the stellar models used for this study, the overshoot 
has been parametrized so that the convective core is extended by 
0.1 H p for stars more massive than 1.7 Mq, 0.05 Hp for stars be¬ 
tween 1.25 and 1.5 Mq and none below. This calibration has been 
made for the rotating models to reproduce the width of Galactic 
open clusters (see Ekstrom et al. 2012). Non rotating models use 
the same calibration for the overshoot. Thus the differences seen 
in Fig. 1 results only from the effect of rotational mixing. 
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Figure 2. Comparison of isochrones in the theoretical HR di¬ 
agram, for age 350 Myr, Z=0.006 and different stellar angular 
velocity. Looking at the MSs, from left to right we have, in units 
ofta cr i t , w=0 (dashed blue), 0.6 (dotted black), 0.8 (double dashed 
black) and 0.9 (full line red). 

the eMSTO can not be due to photometric errors, differen¬ 
tial reddening or unresolved binaries. The authors show that 
if the stars in NGC1856 are chemically homogeneous and the 
eMSTO is due to a prolonged star formation, the star for¬ 
mation must have lasted about 150 Myr (not necessarily in 
a continuous way). 

As discussed in §2 a range of stellar rotation velocities 
can mimic an age spread; here we show which distribution 
of rotational velocities is needed to reproduce the eMSTO 
and of the split of the upper MS. 

In Fig. 3 we show the comparison between data and rotating 
versus non-rotating stellar models. The isochrones, both for 
350 Myr, are taken from the Geneva database and are calcu¬ 
lated from the models by Georgy et al. (2013). The metal¬ 
licity we adopt is Z=0.006 that, the metallicity, available 
on the database, closest to that inferred from spectroscopy 
(Harris & Zaritsky 2009) for the young populations of the 
Large Magellanic Cloud. 

For the non rotating case we extend the isochrones to the 
masses M<1.6Mq by means of the models published by 
Mowlavi et al. (2012). The extension helps to choose a suit¬ 
able reddening and distance modulus, which are labelled 
in Fig. 3. In order to reproduce the split in the upper 
MS, for the selected age, the population of rotating stars 
must be characterized by a rotation rate close to break-up 
(oj ~ 0.9 x u; C rit) (Fig. 3, left panel). Notice that the split 
of MS is less evident in the color mF 43 sw-niFsi 4 w (Fig. 3, 
right panel), and the MSTO has a vertical structure. For 
this reason, UV photometry is crucial to shed light on the 
eMSTO phenomenon in NGC 1856. With this choice of ro¬ 
tation and age we then select from the database the corre¬ 
sponding synthetic photometry with a random viewing angle 
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m F336W m F555W m F438W m F814W 

Figure 3. Two isochrones of 350Myr for Z=0.006, used for the final simulation shown in Fig. 4, are overimposed to the data from Milone 
et al. (2015). The left panel shows the mysssw vs mF336W _m F555Wi and the right panel shows the mF555W vs niF438W _m F8l4W color 
magnitude diagrams. Two different value of angular velocities (uj = 0.9 X cj (:r j T and ui = 0) are plotted. The non rotating isochrone is 
prolonged towards low masses using the models by Mowlavi et al. 2012. This allows to choose the value of distance modulus and reddening, 
labelled in the figure. The absorptions adopted in the different bands are A 336 W /E(B-V) =5.10, A 438 W/E(B-V) =4.18, A 555 YY /E(B-V) 
=3.27, A 814W /E(B-V) =1.86 from Milone et al. 2015. 


distribution, using the gravity darkening model by Espinosa 
Lara & Rieutord (2011) and including the limb darkening 
effect (Claret 2000). The data of the synthetic simulation 
were transformed into observational planes using the model 
atmospheres by Castelli & Kurucz (2003), convolved with 
the HST filter transmission curves. 

In the simulation we take into account the combined pho¬ 
tometry for 30% of binaries, following the procedure out¬ 
lined in Milone et al. (2009). Binaries are included both in 
the rotating and non rotating group. The mass function of 
the secondary stars is randomly extracted from a Salpeter’s 
like mass function, with a lower limit of 0.4 Mq. Notice that 
only the photometric consequences of the presence of such 
binaries are monitored in the simulations (for the evolution¬ 
ary consequences, see the discussion in §4.1), and that the 
choice for the percentage of binaries is chosen based on the 
estimate by Milone et al. (2015). 

As shown in Fig. 4, this simulation describes well the eM- 
STO and the bimodality of the upper MS. 

In other aspects the simulation is not completely satisfac¬ 
tory. In particular it fails to reproduce the data after the 
end of the central H-burning phase. We suggest that this 
may depend both on the parametrization of the overshoot 
in the Geneva models and on the way in which the luminos¬ 
ity and temperature are corrected for the inclination angle 
(see also discussion in §4.1). 

The simulation fails also to reproduce the red clump (panel 
c in Fig. 4), but for a possible solution of this problem see 
§4.1. 


4 DISCUSSION 

4.1 General adequacy of rotating models 

The results shown in the panels of Fig. 4 rely on the Geneva 
database, and should be considered as a first attempt to 
explain the eMSTO and the split upper main sequence in 
NGC 1856 in terms of differences in the rotational velocities 
of stellar subsets. Noticeably, the Geneva database was built 
before the recent discoveries on NGC 1856 and it is remark¬ 
able how the models do such a good job in reproducing its 
split main sequence. 

There are two incomplete issues in the synthetic mod¬ 
els: first, the red side of the eMSTO region is not well re¬ 
produced. It looks like the end of core hydrogen burning is 
not as sharply defined in reality as it occurs in the mod¬ 
els we have used. This looks plausible, as the rotation and 
mixing issue may be subject to several minor second-order 
parameters which can not easily be taken into account in 
modeling. Second, the splitting of the main sequence is best 
reproduced if we assign a very fast rotation (w=0.9 w C rit) to 
the numerous sample (67% of stars) defining the redder MS 
and the upper turnoff. This may depend on the way gravity 
and limb darkening are introduced in the models. This is 
the first adequate representation of the splitting, so it is a 
powerful hint that rotation is the reason for the CMD fea¬ 
tures, even if we should not take the very high rotation at 
face value. 

Alternatives to the rotational interpretation can reside in 
different metallicity, or different specific elemental abun¬ 
dances for the two sides of the main sequence, but further 
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Figure 4. rn F555W vs rnj- 33 g\y- mp 5 55 yy color magnitude diagram of cluster NGC1856 and its Hess diagram. Data (grey open circle) 
from Milone et al. (2015) are compared with the simulations (filled triangles) described in detail in the text, at the luminosity of the 
turnoff (panel b); the red clump (c); and the upper part of MS (d). For this last case, the histogram of the A (mp 33 q yy-rnI 555 VV) 
distribution for stars in six magnitude intervals is reported in panel e, after having verticalized the upper part of the main sequence 
delimited by the gray boundaries in panel b). 


modelling is required, and whether such differences exists 
should be checked by future spectroscopic observations. 

The final problem is that the clump reproduction is 
quite modest. While the rotating simulation reproduces the 
luminosity location of the clump, there should be a second 
clump (less populated) about 1 magnitude dimmer, but the 
data do not show it. The clump evolution is a very delicate 
issue, the ratio of blue to red stars depending on the de¬ 
tailed treatment of convection and mixing. Consideration of 
non instantaneous mixing produces favors a long evolution 
in the blue side of the clump (e.g. Ventura et al. 2005), so it 
could describe better the clump here. Here again, notice that 
age differences produce a similar problem for the clump(s) 
location. A possible solution of this problem occurs natu¬ 
rally if the slow-rotating stars are all binaries, as discussed 
in §. 4.2. Observational data and theoretical considerations 
show that orbital periods between 4 and 500 days can be 
effective in producing slow stellar rotation (Abt & Boon- 
yarak 2004). As we are dealing with evolving masses ~3 Mq, 
orbital periods between ~50 and 500 days, according to Ke¬ 
pler’s law, would produce Roche lobe interaction and mass 
transfer from the primary star while it is evolving along the 


RGB. Shorter periods produce mass transfer at an earlier 
phase, during the crossing of the Hertzsprung gap. Thus the 
evolving star would not reach helium ignition as a single 
star would, and the low luminosity blue hook could result 
depopulated. In the end, the lack or sparseness of a low lu¬ 
minosity clump might lend further support to the hypothesis 
of a rotational origin for the splitting of the turnoff. 

A consequence of this interpretation is that the role of 
binaries in the CMD should be reconsidered. The slowly ro¬ 
tating group should be made of 100% binaries with periods 
between 4 and 500 days, while the rapidly rotating group 
would include both single stars and longer period binaries. 
The binaries clearly seen on the right of the observed main 
sequence(s) could be mostly belonging to the non rotating 
population. A full discussion of this problem goes beyond 
the purpose of this work, and is delayed to a following in¬ 
vestigation. 


4.2 Why two populations? 

Our analysis has shown that in order to reproduce the ob¬ 
served features of the CMD of NGC 1856 it is necessary to 
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have a large fraction of cluster stars (2/3) retaining the large 
rotational velocity typical of the B and A stars (see e.g. Zorec 
& Royer 2012, for A-type stars), while the remaining 1/3 of 
the stars must have slowed down early enough to prevent ro¬ 
tational mixing, and have evolved as standard non-rotating 
or slowly rotating models. 

The key question to address therefore concerns the 
mechanism and processes responsible for slowing down a 
large fraction of stars during the H-core burning evolution. 

Zorec & Royer (2012) have gathered and analyzed the 
rotational velocities of fields A-type main-sequence stars 
(types B6 to F2). They found that stars below 2.5Mq show 
an unimodal velocity distribution displaying a small accel¬ 
eration as the stars evolve along the MS. On the contrary, 
stars between 2.4 and 3.85M 0 display a clear bimodal dis¬ 
tribution with a slow velocity component (up to 20% for the 
2.8M 0 bin). They found that the two peaks are clearly sepa¬ 
rated during the MS evolution but these peaks tend to bend 
near the turn-off. This dichotomy in velocity persists if close 
binary and chemically peculiar stars are included. The origin 
of the slow component is unknown but the authors make two 
hypotheses: (1) some stars have already lost all their angu¬ 
lar momentum during the pre-MS due to magnetic braking 
or (2) in a fraction of undetected binaries, tidal interactions 
have slowed down the stars. It should also be noted that 
a similar bimodality in the velocity distribution have been 
detected among early B-type stars by Dufton et al. (2013). 

Abt & Boonyarak (2004), in a study of the velocities 
of A and B-type binaries, found that close binaries (period 
smaller than 4 days) are all synchronized and circularized. 
Binaries with period longer than 500 days display the same 
rotational properties as single stars. For binaries with period 
between 4 and 500 days, rotational velocity is significantly 
smaller than for single stars, with about one-third to two- 
third of their angular momentum being lost, presumably, by 
tidal interactions. 

On the theoretical side, tides in binaries will be respon¬ 
sible to the synchronization and circularization of the orbit 
(Zahn 1977). For low-mass stars with a convective envelope, 
the main mechanism is the viscous dissipation of the ki¬ 
netic energy acting like an equilibrium tide. On the other 
hand, in more massive stars with a convective core, some low 
frequency modes of oscillation can be excited by the peri¬ 
odic tidal potential. The response to this excitation is called 
the dynamical tide (see e.g. Kopal 1968). With a simplified 
model, assuming an uniform rotation, Zahn (1977) obtains 
the following expression for the synchronization time: 


1 


tsync 


5 x 2 5/3 



1/2 

q 2 (l + <l ) 5/6 




(1) 


where q is the mass ratio of the binary, f ?2 is the parame¬ 
ter measuring the coupling between the tidal potential and 
the gravity mode which depends on the convective core mass 
(~ 5 x 10 -8 for a 3M 0 ). Zahn (2008) shows that the synchro¬ 
nization time scales with (7 /MR 2 )/E 2 R 7 , which increases 
during the main sequence. Thus, we can conclude that most 
of the tidal interaction will take place around the ZAMS, 
so that binaries dispersion, if it occurs during the main 
sequence, will take place in systems already slowed-down. 
Moreover, early-type binaries have their transition period 


(for circularization) observed through OGLE and MACHO 
survey well in agreement with the dynamical tides theory 
(North & Zahn 2003). However an important fraction of bi¬ 
naries are also circularized even if the semi-major axis is 
larger. This may indicate that another more efficient mech¬ 
anism is at play for them. One possibility will be that these 
binaries have undergone several episodes of resonance lock¬ 
ing (see Witte & Savonije 1999b,a, for the the effect on a 
1OM 0 primary) 

In summary, there are several ways to slow down a frac¬ 
tion of the cluster stars (see also the extended discussion in 
Dufton et al. 2013). If we look into more detail into the 
binary hypothesis, we should expect that most of the non¬ 
rotating stars hide a binary companion. In fact, binaries of 
periods in the range 4-500 days would be hard binaries in 
the core of this cluster, whose central dispersion velocity is 
~2.5 km/s (McLaughlin & van der Marel 2006). Only bi¬ 
naries with period longer than ~ 10 3 yr would be ionized in 
the central regions. Notice that the cluster seems to host 
~30% (Milone et al. 2015) of binaries, so this hypothesis is 
not unreasonable. Further investigations may provide fur¬ 
ther insight into this problem. 


4.3 The case of multiple populations in old 
globular clusters 

We conclude our discussion with some comments and cau¬ 
tionary remarks about the connection between eMSTO clus¬ 
ters and multiple-population old GCs. 

Evidence of multiple populations in old GCs comes from 
an extensive number of photometric and spectroscopic stud¬ 
ies showing that stars in old GCs are characterized by a 
spread in the abundances in helium and light elements (such 
as Na, O, Mg, Al). Photometric studies have also shown the 
widespread presence of discrete groups with distinct photo¬ 
metric properties corresponding to different chemical prop¬ 
erties (see e.g. Piotto et al. 2015). These are the key finger¬ 
prints of multiple populations in old GCs and provide the 
main constraints for any theoretical effort aimed at modeling 
the formation history of these systems and shedding light on 
the source of processed gas composing the second-generation 
populations observed in globular clusters. 

While the study of young and intermediate-age clusters 
(hereafter YICs) might provide some clues to the answers 
to the many questions raised by the discovery of multiple 
populations in old GCs it is important to exercise much 
caution in connecting YICs and old GCs for this purpose and 
in drawing conclusions about different theoretical scenarios 
for the formation of old GCs from observations of YICs. 

Some specific issues to consider are the following. 

1) It is not clear whether the YICs studied so far ac¬ 
tually host multiple populations at all. As mentioned in the 
Introduction, for example, a spectroscopic study of one of 
the eMSTO clusters (NGC 1806) does not reveal the typ¬ 
ical chemical fingerprints of multiple populations found in 
old GCs. Photometric signatures which might be ascribed 
to extended star formation history in YICs, such as those 
discussed in this and other works, have possible alternative 
explanations unrelated to the presence of multiple stellar 
populations with different chemical properties. It is impor¬ 
tant to realize that spectroscopic characterization of YIC 
aimed at shedding light on the presence (or lack of) the 
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chemical patterns typical of multiple populations is desirable 
to make a solid connection between YIC and old multiple- 
populations GCs even if strong evidence of an age spread in 
YIC were to be established. Although, as remarked in §.l, 
many additional spectroscopic studies are necessary to clar¬ 
ify this issue, it is possible that the study of these clusters 
might be irrelevant to shed light on the formation history of 
multiple-population old GCs. 

2) While eMSTO, as also argued in this paper, might 
indeed be explained without relying to extended periods of 
star formation (and eMSTO clusters might not host chemi¬ 
cally different multiple populations), studies ruling out star 
formation episodes extending for a few hundred Myr (see, 
e.g. Bastian & Silva-Villa 2013; Bastian et al. 2013) have 
been presented as evidence against the AGB scenario of 
multiple-population formation. We point out here that in all 
the studies (see e.g. D’Ercole et al. 2008, 2010, 2012) aimed 
at modeling the chemical properties of multiple-population 
old GCs using AGB ejecta the epoch of second-generation 
formation is limited to a time interval between about 30 and 
100 Myr with a significant fraction of second generation stars 
forming in the early phase of this period. While in a few clus¬ 
ters second-generation star formation might have extended 
beyond (see e.g. D’Antona et al. 2011) the current observa¬ 
tional constraints from spectroscopic studies - in particular 
the lack of C+N+O increase in the second generation stars of 
most clusters— do not, in general, allow such an extended 
period of second-generation star formation. In addition to 
urging caution in connecting YICs and old GCs, it is impor¬ 
tant to realize that lack of evidence of such an extended age 
spread can not be used as evidence against AGB models. 

3) Finally we emphasize the importance of properly 
considering the differences between the structural proper¬ 
ties of YIC and old GCs. The current properties of old 
GCs have been affected by early and long-term evolution¬ 
ary processes (e.g. gas expulsion, mass loss due to stellar 
evolution, two-body relaxation, tidal interactions). Recon¬ 
structing their structure and mass at the time of second- 
generation formation in order to identify the threshold in 
the structural properties allowing the formation of multiple 
populations is not a trivial task. Identifying young clusters 
with present masses similar to the current masses of old 
multiple-population clusters to test multiple-population for¬ 
mation scenarios, without considering the complications of 
properly reconstructing the evolutionary history of old GCs, 
may easily lead to erroneous conclusions. For example the 
current mass of NGC 1856 is about 10 5 Mq (McLaughlin 
& van der Marel 2005), similar to the current mass of some 
old GCs hosting multiple populations. Its structural prop¬ 
erties are such the central escape speed is about 10 km/s 
(McLaughlin & van der Marel 2005). Detailed models would 
be necessary to understand whether initial/earlier properties 
allowing the formation of a significant fraction of second- 
generation stars are among those evolving into NGC 1856’s 
current properties —see Correnti et al. (2015) for an attempt 
to address this issue; see also Goudfrooij et al. (2014). How¬ 
ever the current structural properties, if sufficiently similar 
to those of the cluster at earlier stages of its life, would in¬ 
dicate that this cluster is not capable to retain a significant 
fraction of the ejecta of AGB or other possible sources of 
processed gas, lend support to the interpretation of the eM¬ 
STO in terms of different rotational velocities, and provide 


an example of the caution necessary in using these stellar 
systems to draw any conclusion about multiple population 
formation models. 
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